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Abstract 

We discuss inflaton decay in supergravity, taking account of the gravitational effects. It is shown 
that, if the inflaton has a nonzero vacuum expectation value, it generically couples to any mat- 
ter fields that appear in the superpotential at the tree level, and to any gauge sectors through 
anomalies in the supergravity. Through these processes, the inflaton generically decays into the 
supersymmetry breaking sector, producing many gravitinos. The inflaton also directly decays into 
a pair of the gravitinos. We derive constraints on both inflation models and supersymmetry break- 
ing scenarios for avoiding overproduction of the gravitinos. Furthermore, the inflaton naturally 
decays into the visible sector via the top Yukawa coupling and SU(3)c gauge interactions. 
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I. INTRODUCTION 

Inflation [1] provides a simple solution to a number of serious shortcomings in the big 
bang cosmology such as the horizon and flatness problems. Above all, it can account for the 
origin of density fluctuations necessary to form the rich structure of our universe. In fact, 
the Standard slow-roll inflation predicts almost scale-invariant power spectrum, which fits 
the recent cosmic microwave background (CMB) data [2J quite well. 

It is now recognized that the universe underwent an inflationary epoch at an early stage. 
During the inflation, the universe is dominated by the potential energy of the inflaton, and 
experiences exponential expansion After inflation ends, the inflaton field releases 

its energy into a thermal plasma by the decay, and the universe is reheated. Since all the 
particles including photons and baryons in the present universe are ultimately originated 
from the inflaton decay, it is of great importance to reveal how the reheating proceeds. 

So far however, the reheating process has not been fully investigated. One often sim- 
plifies the whole reheating processes, and expresses them in terms of a single parameter, 
the reheating temperature. That is, the inflaton is assumed to have some ad hoc interac- 
tions with lighter degrees of freedom, i. e., the Standard model (SM) particles in most cases, 
while possible productions of the hidden fields and/or gravitinos are neglected without def- 
inite grounds. However, many cosmological phenomena, e. g., baryogenesis, and production 
of dark matter and unwanted relics, crucially depend on the details of the reheating. Al- 
though the reheating temperature is certainly an important characteristic parameter, such 
simplification is too crude to truly describe cosmological scenarios. 

Recently there has been much progress concerning the decays of scalar fields such as 
moduli 0, 0, (]] and inflaton lloL fll| in a framework of the local supersymmetry 

(SUSY), i. e., the supergravity (SUGRA). The supersymmetric extension is one of the most 
promising candidates for the theory beyond SM. If SUSY exists at the TeV scale, the inflaton 
dynamics is quite likely described in SUGRA. In addition, since the existence of a flat 
direction is mediocre in SUSY models, one can find extremely flat potentials appropriate 
for the slow-roll inflation. Throughout this paper we consider inflation models in SUGRA. 
We have investigated the reheating of the universe in this framework, and found that the 
gravitinos are generically produced from the inflaton decay in most inflation models. In 
particular, Ref. [7J has first pointed out that the inflaton can directly decay into a pair of 
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the gravitinos. Moreover, incorporating the gravitational effects, Refs. [9|, [llj have shown 
that the inflaton generically decays into the SUSY breaking sector, which produces the 
gravitinos (in)directly. The gravitino production rates due to these processes depend on the 
inflaton parameters as well as the detailed structure of the SUSY breaking sector. Such 
gravitino production clearly goes beyond the simplification of the reheating that has been 
adopted so far, and interestingly enough, it provides severe constraints on inflation models 
as well as the SUSY breaking scenarios. These constraints, together with the future collider 
experiments and observations on CMB, should become an important guide to understand 
the high energy physics and the early universe. The purpose of the present paper is to 
provide a global picture of the inflaton decay processes in SUGRA, paying special attention 
to the gravitino production. In particular, we explain which decay processes become most 
important under which circumstances. Not only do we summarize the decay processes found 
so far but we also give complete results on the spontaneous decay and the anomaly-induced 
decay processes, including the higher dimensional terms in the Kahler potential and the 
Kahler and sigma-model anomalies. 

The organization of the paper is as follows. In Sec. [TT] we review the gravitino pair 
production at the inflaton decay. Then we discuss the spontaneous decay at the tree level in 
Sec. IIHI In Sec. |IV]we consider the anomaly-induced decay of the inflaton, which proceeds via 
the anomalies in SUGRA. We provide some results on the decay rates, by way of illustration, 
for the minimal and sequestered Kahler potentials there. In Sec. IVj we study cosmological 
implications of the decay processes explained in the preceding sections, particularly focusing 
on the constraints on the inflation models. The last section is devoted to conclusion. 



II. DECAY INTO A PAIR OF GRAVITINOS 



Once the inflaton field obtains a finite vacuum expectation value (VEV), it necessarily 
decays into the gravitinos. In this section, we briefly discuss the production of a pair of the 
gravitinos, and provide the partial decay rate. The process we consider is a perturbative 
decay, and the gravitinos are produced directly from the inflaton. The gravitino production 
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is represented by the following interactions in the SUGRA Lagrangian [U 

~e G ' 2 (GiP + Grt**) [^ m o mn ^ n + $ m v mn $n] , (1) 

where a mn = \{a m a n — a n a m ), and we have chosen the unitary gauge in the Einstein frame. 
The sum over the indices is understood unless otherwise stated. We have also adopted the 
Planck unit Mp = 1 (Mp = 2.4 x 10 18 GeV) here and in what follows unless it is written 
explicitly. The 2-spinor, ip m (or ^3/2), represents the gravitino, while 0* collectively denotes 
an arbitrary scalar field including the inflaton 0. Then the decay rate of the inflaton into a 
pair of the gravitinos, ~^ 2-0 3 / 2 , is evaluated as [4| 

F (grav) _ \ G 4>? m \ / 2 x 

" 288tt m 2 /2 M 2 ' K ' 

where ms/2 = e G ^ 2 and are the masses of the gravitino and the inflaton, respectively. We 
readily find that the decay amplitude is inversely proportional to 777.3/2- This is a result of an 
enhancement (oc m^? 2 ) due to the longitudinal mode of the gravitino, ip m (k) oc k m /m 3 / 2 ~ 
777^/777,3/2, which is partially compensated by the chirality suppression of the amplitude 
(oc 7773/2). 

The decay amplitude crucially depends on G^, which is a derivative of the generalized 
Kahler potential, G = K + ln \ W\ 2 , with respect to the inflaton field 0. It is related to an 
F-term of the inflaton supermultiplet through the equation of motion, F 1 = —e G ^ 2 g^*Gj*. 
In order to evaluate G^, we need to incorporate the SUSY breaking field, z, into our analysis. 
This is because of the following reason. The decay is treated in the mass-eigenstate basis. 
In this basis, <fi generally mbces with z due to the SUGRA effects, unless the inflaton is 
protected by some symmetries which are preserved at the vacuum. We take G z = 0(1) to 
have the vanishing cosmological constant. Then, G z can contribute to G^ effectively via 
mixings between and z, which enhances gravitino production rate from inflaton decay. 
That is, the inflaton first oscillates into z^ , which then decays into a pair of the gravitinos: 

# Z« -> 2^3/2- 



Due to the Kahler invariance, the generalized Kahler potential G is more convenient and transparent 
than using the Kahler potential K and the superpotential W. Since these two frames are related by the 
Weyl transformation, any physical amplitudes are equivalent at the tree level. 
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The mixing angle depends on the mass spectrum of and z. The direct pair-gravitino 
production is effective especially for <C m, , Such a large m z is often realized in the 



dynamical SUSY breaking (DSB) scenario 13j. In this case, there is a soft mass term, 
K ~ |z| 4 /A 2 (A is the dynamical scale), and a scalar mass of z can be larger than m^, 



depending on inflation models. Then G& is given by 



6] 



\G^ ~ V3g</>z* 



2 

+ 



V3(V G 2 



m 3/2 2 



(3) 



where we have neglected interference terms and higher dimensional operators in the Kahler 
potential. Here g^* = d ^^j and V^Gj = Gij — T^Gk with T l - k = g u * gje*k- Note that the 
first term is from the mbdng in the kinetic terms, while the SUGRA effects contribute to 
the second one. Tirus even if there is no direct coupling between the inflaton and SUSY- 
breaking sectors in the global SUSY limit, the inflaton decays into a pair of the gravitinos 
for <C m z . 

Since each term of (j3J) is expected to depend on linearly, it is convenient to express the 
mixings as 

\V^G Z \ = c(<f>), 
\g^z*\ = c (cp). (4) 

In SUGRA, c is estimated to be 0(1) for a generic Kahler potential by using G z = 0(1), 
while c depends on details of the SUSY breaking sector such as the VEV (z), e. g. c = (z) 
for 8K = \<f)\ 2 \z\ 2 . Then, if c is suppressed as in case of the minimal Kahler potential (i. e. 
g<j>z* — 0), the gravitino pair production rate is 



96tt \M p Ml' 



On the other hand, if the kinetic mbdng is large, the rate is much enhanced as 



c 2 /wy ml 

96tt \Mp) ml /2 M'p' 



(6) 



Such large gravitino production rates are cosmologically disastrous, which will be discussed 
in Sec. EVJ 

For high-scale inflation models with 3> m z , the pair-gravitino production rate depends 
on the detailed structure of the SUSY breaking models. If the SUSY breaking field is 



singlet and elementary above A b c , the inflaton still directly decays into a pair of the 
gravitinos. In this case the relevant contribution to G$ comes from higher dimensional 
terms, K ~ (K,/2)\(f)\ 2 zz + h. c. d . Then the gravitino production rate is given by ([S]) with c 
replaced with k (see Ref. jfj, 7] for details). On the other hand, if the SUSY breaking field is 
composed of other fields and if the dynamical scale A is below m^, the direct production of the 
gravitinos becomes suppressed. Instead, as discussed in the following sections, gravitational 
effects force the inflaton to decay into the SUSY breaking sector. 

Finally let us make a comment. In addition to the pair-gravitino production, the gravitino 
may be singly produced at the decay. This is the case when the inflaton mass comes from 
the soft SUSY breaking terms. Then the rate becomes as large as that given by with 
c = 0(1). However, if the inflaton mass is provided by a SUSY-invariant mass term (as in 
most inflation models), such a single-gravitino production is negligible. 

III. SPONTANEOUS DECAY 

In this section we review the spontaneous decay of the inflaton, 0, at the tree level. 
If the reheating is induced by the inflaton decay through non-renormalizable interactions, 
the reheating temperature can be low enough to satisfy the constraints from gravitinos 



produced by thermal scatterings [18|, ll9|]. Since the interactions are then quite weak, the 
SUGRA effects may play an important role. Indeed, it has been recently pointed out that 
the SUGRA effects induce the inflaton decay The relevant channels of the inflaton decay 
contains the 2- and 3-body final states e . 

For the matter-fermion production, the relevant interactions are provided in the Einstein 



frame as 



12| 



b Such a singlet SUSY breaking field is necessary for the gauginos to have a sizable mass in the gravity- 

mediated SUSY breaking scenario 0, 0|- See also Ref. [16] for a retrofitted gravity-mediation. 
c In the DSB scenarios, the Polonyi problem was once solved since the z field decays much before BBN 

due to a large soft scalar mass of z HIU. However, it has been recently found that the presence of such 



a field still puts a severe bound on the inflation scale 17 1. 
d The contribution from this operator is suppressed when is smaller than m z 
c Although (j) may decay into 4-scalars, it is suppressed by the phase space and can be neglected. 
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-I e ^(p^. W ) x Y+h.c, (7) 

where ViDjW = Wij + K {j W + K, L DjW + KjDiW - K { KjW - T%D k W with D t W = 
Wi + KjVV '. On the other hand, the matter-scalar production is represented by the kinetic 
term and the scalar potential; 

e' 1 C = -g^d^d^ 

'g^^DiW^DjWy - 3\W\ 2 ] . (8) 



-e K 



In this paper, we adopt notation that, when the inflaton is explicitly shown in expressions, 
4> l and \ l represent only the matter fields. (We also use Q to denote the matter fields.) Oth- 
erwise, as mentioned before, (f) 1 collectively denotes an arbitrary field including the inflaton 
0. 

First let us consider decay processes induced by higher dimensional operators. The rel- 
evant one arises from such terms in the Kahler potential that holomorphically depend on 
the matter fields, Q, Q', i. e. 5K ~ \<p\ 2 QQ' + h. c. f . The presence of such operators strongly 
depends on symmetries of the visible/hidden sectors. In Eq. ((7j) the term including Tfj in 
the fermion mass is given by 

C = V/V^*^0xV + h.c., (9) 

which induces the inflaton decay into the two fermions (<f) — > x*X"0- Note here that 
e K l 2 g^*W^ = is the inflaton mass. 

On the other hand, the decay into the two scalars (0 — > i 0- J ) arises from the kinetic term 
of the matter scalars, 

£ = ^^(9 2 0)0*>* J +h.c. (10) 

Using the equation of motion, d 2 (p = m^cp, one finds that the decay rates satisfy F(<p — > 
~ T(0 — > The total rate then becomes 



r (2-body;hol) _ p^ ^ -i -j j + p 



8ir \ mi 



The decay process from this operator is obtained also in the global SUSY models. 
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where Mq is a mass of the final state particles. Here note that i and j are fixed and the 
sum is not taken over these variables in the last expression. 

Next let us discuss the case of Tfj = 0, which is due to some symmetries imposed on 
the i- and j-matter fields. The 2-body decay then becomes suppressed by the mass of the 
final-state particles, Mq. From the Lagrangian (17j) , the effective interaction is given by 

C = -\e K ' 2 (K+Wtj + - 2Y%W jk ) X ¥ + h.c., (12) 

where we have assumed that |I?0W| <C \W\. Here and in the followings, we assume that the 
matter fields are charged under some symmetries for simplicity. Then we can set K il Wi -C 1 
for the matter fields. It should be noticed that the second term in the bracket is necessary 
to ensure the Kahler invariance. For instance, if we apply the Kahler transformation, K — ► 
K — (K^)(j) — [K (j,*) (j)*, the first term vanishes and the second term compensates it. This 
becomes clear if we write the interactions in terms of G. The effective Lagrangian f[T2"j) is 
represented as 

C = -\e G l\G^ - 2r* i G i *)foV + h.c, (13) 

which is obviously invariant under the Kahler transformation. Note that rL in Eq. f|T2|) 
is different from Vfj in Eq. (J9]). The coefficient, ~ K^*, can be nonzero easily. For 
instance, 5K ~ |0| 2 |Q| 2 leads to T®q ~ (0), which is nonzero as long as (<p) ^ 0. 

On the other hand, the 2-scalar production consists of the two channels; — > cj) 1 ^ and 
(f) — > The former comes from the scalar potential; 

C = - V (KfWij + W*, - 2T%W jk y g^W^ + h.c. (14) 

We can easily check that this provides the same decay rate as that of the fermion final state 
induced by (jT2l . i. e., T(0 — > x l x J ) = T(0 — >■ l J ). Also the kinetic term of the scalar 
fields gives another decay channel, <fi — > However its amplitude is proportional to the 

scalar mass squared of the final state, noting [0|<5| 2 ] D = 4>{d 2 Q*)Q + ■ ■ -. Thus the process 
becomes dominant only when the scalar fields has a quite large soft scalar mass. 

To summarize, the total decay rate of the 2-body final state from the interactions (fT2l) 
and flUD is 

r (2-bod y ) s r^^^+r^^^) 

m A l -^Jr) . ( 15 ) 



'% 3 



8n \ m 




FIG. 1: The decay of the inflaton into the three-body final states; the decay with the four-point 
vertex, and with the fermion and scalar exchanges, from left to right. 



where C\f = e K \K 4 ,W ij + - 2r^I47 jfc | 2 with fixed i and j (the sum is taken only over 
k). If the particles in the final state have a SUSY mass, W = M Q QQ' ) C\f is proportional 
to Mq. If the two particles in the final state are identical to each other, W = \MqQQ (e. g. 
the right-handed neutrino N with a Majorana mass M/v), the decay rate becomes half of 

(USD s. 

Next we consider the decay with 3-body final states. The decay processes through the 
dimension five operators are — > <p* l X^X k an d <P ~ > (j) 1 $ (j) k . The former process is composed 
of the three diagrams in Fig. UJ In addition to the spontancous decay process pointed out 
in Q] (the left diagram), the higher dimensional terms in the Kahler potential contribute 
to the decay rate (middle and right). Evaluating these diagrams, we obtain the effective 
interaetions as 

C = - l -e Kl2 (K^W ijk + W^ k - 3I%W jM ) + h.c. (16) 

On the other hand, the interaetions representing the decay into 3 scalars, <f> — > (p' 1 'cfp (f) k , are 
obtained by expanding the scalar potential as 

C = - X -e K (K^W ijk + W 4>ljk -3T^W jke y g^'W^^^^ + h.c. (17) 

One can write down these interaetions in terms of the Kahler invariant funetion, G, by 
replacing K^Wijk + W^jk — > G^jk and Wju — > Gj k g, respectively. We find that the decay 



s The spontaneous decay into the right-handed (s)neutrinos make the non-thermal leptogenesis seenario 
quite attractive 



20] 
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rate into 3 scalars is same as that into 1 scalar + 2 fermions, i. e., T(0 — > 0*0- J fc ) ~ T(0 — > 
< f ) * l X^X k ) + r(0 — > (j>*^X k X l ) + r(0 - ► 0* fc X*X"O> f° r fixed i, j and /c. Summing these decay 
rates, the total 3-body decay rate is given by 

r (3-bod y ) = r(0 ^ 3 scalars) + T(0 -> 1 scalar + 2 fermions) 

r (3) 

where = e K \K < pWijk + W^jh — S^jWikif' with fixed i, j and (sum over f). Here we 
have neglected the masses of the fmal-state particles. 

Finally we discuss the inflaton decay into the gauge bosons and gauginos. At the tree 
level, it is effective only when the gauge kinetic function depends on the inflaton field h . 
Actually, we obtain the total rate of the decay into the gauge sector as 4] 

r (gaugetree) „ (lg) 
47T 

from C — k J d 2 9 (pW a W a , where W a is a field strength of the gauge supermultiplet, N g 
is a number of the generators of the gauge symmetry, and we have assumed the canonical 
normalization for the inflaton and gauge multiplet. In f|T9|) . half of the decay rate comes 
from the gauge boson production and the other half is from the gaugino production. 

Except for such direct couplings, no sizable interactions are found at the tree level be- 
tween the inflaton and the gauge fields in the SUGRA Lagrangian [l^j]. The feature can be 
understood by using the gravity supermultiplet. The multiplet is minimally composed of 
the followings; 

h mn , <, b m , M, (20) 

which represent the graviton, the gravitino, and the vector and scalar auxiliary fields that 
correspond to the U(1)r and conformal symmetries of the superconformal transformation, 
respectively. Even in the absence of the direct couplings, the gravity multiplet can connect 
the inflaton field to the visible/hidden sectors. Actually, the auxiliary fields, b m and M, 
depend on the inflaton field as well as the visible/hidden fields due to the equation of motion, 
and the longitudinal component of the graviton, h, is related to the inflaton through the 
Lagrangian term, C = —^e~ K / 3 7Z, in the conformal frame \ 



h The coupling may be induced by the mixing of the inflaton with other fields such as the SUSY breaking 
field 



1 Even in the non-SUSY models, the latter contribution can arise [21 1. 
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The relevant terms involving b m in the SUGRA Lagrangian 12] are given by 

C aux = h m b m - ^(KM* - A^«9 m 0*> m + l -g i3 *x^mX l b m - \\a m Xb m . (21) 

Solving the equation of motion for b m , one can see that b m depends linear ly on with a 
coefficient K^, and that it also includes the gaugino current, A<r m A. The decay into a pair of 
the gauginos is thus suppressed by the gaugino mass because the processes requires a chirality 
flip. In other words, noting that the inflaton contributes to the longitudinal component of 
b m , the U(1)r charges of the final state should be nonzero for the decay to proceed due to 
the U(1)r current conservation. Thus the gaugino mass appears in the amplitude. 



221 is 



Next we focus on h and M. The superconformal calculus formulation of SUGRA 
convenient to understand the decays mediated by these fields. In fact, their contributions 
can be taken into account by incorporating the chiral compensator field into the Lagrangian. 
The F-term of the compensator corresponds to M by using the equation of motion, and M 
includes a linear term with respect to 0, whose coefficient is proportional to K^. Further, 
since the compensator has a Weyl charge, its scalar component depends on K after the 
Weyl transformation to canonicalize the gravity sector, i. e. from the conformal frame into 
the Einstein one. Then <fi arises linearly in the scalar component when is non-zero. Thus 
the operators induced by h and M are represented by multiplying the compensator field. 
It is, however, known that the compensator does not physically couple to the gauge sector 
because it is conformal. Consequently, the decays into the gauge sector are suppressed at 
the tree level. 

Before closing, it is interesting to note that these features are broken at the quantum 
level. That is, the inflaton can decay into the gauge sector via anomalies. We will discuss 
this mechanism in the next section. 



IV. ANOMALY-INDUCED DECAY 

At the classical level, the spontaneous decay of the inflaton into the gauge sector is 
suppressed, since the gauge sector is conformal as discussed in the previous section. The 
quantum corrections, however, violate the conformal invariance, and so, the inflaton decay 
into the gauge sector may arise at the quantum level. Taking account of the SUGRA 
effects, the super- Weyl-Kahler (SW-Kahler) symmetry and the sigma-model isometry are 
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anomalous at the quantum level. Not only do these anomalies mediate the SUSY-breaking 
effects to the visible sector , but they also enable the inflaton field to couple to the gauge 
supermultiplets [111 ]. 

In the superfield description, the 1PI effective action includes the non-local terms corre- 



sponding to the anomalies 
AC = 



24 



m-, 



(16vr) s 



d 2 6 W a W, 



D 2 



A(T R - ZT G )R) 



O 2 

-^-D 2 K + ^D 2 log det K\ 



Tr 
3 



T 3 
d R 



R 



+ h.c. 



(22) 



at the leading order of 1/Mp in the conformal frame. Here D is a covariant derivative of 
the supersymmetry, and g is a gauge coupling constant. The coefficients, Tq and Tr, are 
the Dynkin index of the adjoint representation and matter fields in the representation R of 
dimension d R , which are normalized to N for SU(iV) and 1/2 for its fundamentals. A sum 
over the matter fields is understood. Also K\' R denotes the Kahler metric restricted to the 



representation R. The first term in the bracket of Eq. ( |22l corresponds to the SW anomaly, 

and it is not invariant under the S W transformation. In fact, the superspace curvature R 

I 1 

changes under the SW transformation as [12]; 



5R 



-2(2S 



E)R--D 2 E, 

' 4 



(23) 



where a chiral superfield S is defined so as to rescale the vielbein, SEj^ = (S + S)i?^ f , and 
the last term induces a shift of R. On the other hand, the second and third terms in Eq. (1221) 
arise from the Kahler and sigma-model anomalies, respectively. 



i* 



In the conformal frame, R) is expanded as R) = — |[M* + 6 2 (—lZ/2 + id m b m )} + 
where • • • is irrelevant for the decay. In addition to the auxiliary fields, the Ricci scalar, 
1Z depends on the inflaton field through the kinetic term, C = —^e~ K ^ 3 lZ, which induces 
the mbdng of the inflaton with the longitudinal mode of the graviton. To simplify the 
calculation, let us go to the Einstein frame where the gravity is canonically normalized. To 
this end, we perform the SW transformation with = (p e + \/26xe + d 2 F e defined by 26] 

1 



—K, X e = ~KiX l 
12 o 



K { P - -K ijX l X 3 ■ 



Then the anomaly- induced term becomes 26 ] 



2 

AC E = A£ + ^(T r - 3T g ) J d 2 6 T, E W a W a + h.c, 



(24) 



(25) 
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where the fields in AC are simply replaced by those defined in the Einstein frame K 

Expanding the superfields in terms of the components, one obtains interaction terms of 
the inflaton field to the gauge bosons/gauginos k ; 

C = -^X G (j>(F mn F mn - iF mn F mn ) - ^X G m^*XX + h.c, 

X G = (Tq - T R )K^ + — ^(logdet AT|';),^ (26) 

<1r 

where F mn is a field strength of the gauge field and F mn = e mnkl F^ij1. Here we have also 
used the equations of motion for the auxiliary fields in the Einstein frame; 

b m = ^(KM* - K t ,d m <P* 1 ) + • • • , F i = -e K/2 g^(W, + KjW)*. (27) 

It is noticed that M* = —3e K ^ 2 W* does not induce the decay because of \W^\ ~ rn 3 / 2 (0) 
for the inflaton, (f). The total decay rate from (126]) becomes 



p(anomaly) „ ^^2^ (2g) 

where a is a fine structure constant of the gauge group. Note that half of the decay rate 
comes from the decay into the two gauge bosons, while the other half from that into the 
gaugino pair. 

Let us compare the rate of the anomaly-induced decay (|28|) with that of the spontaneous 
decay at the tree level ffTol) and (TTgj) . We find that all these rates are proportional to |-fQ>| 2 . 
It means that, if the Kahler potential of the inflaton is canonical, the VEV of the inflaton 
field is necessary for the decay to proceed by the SUGRA effects. In contrast, the phase 
space and coupling constants depends on each process. The decay rate into the 2-body final 
state (|T5l) is suppressed by the mass squared, Mg/Mp <C 1. While the rate of the 3-body 
final state f|T8|) is suppressed by the phase space compared to (fT5|) . Instead, the coupling 
constant is given by the Yukawa coupling, W^. Compared to these tree- level processes, the 
anomaly-induced decay takes place at the one-loop level. However, since the final state is 2 
body, i. e. a pair of gauge bosons and gauginos, its rate is not negligible compared to those 
of the spontaneous decays at the tree level. 



j A factor in front of Tr is diffcrcnt from the result in 26( because here K in AC is not shiftcd 



k This result is also obtained at the component level by the Weyl rescaling, — ► e~ 2cr e° n , from the 
conformal frame to the Einstein frame. Then the TZ and M shift as 81Z = 12 d 2 a and SM = —KiF 1 with 
a = K /12, while b m remains unchanged. 
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Let us comment on a mass spectrum of the matters in the visible/hidden sectors. In 
this section, we have discussed anomalies that connects the inflaton with the gauge sector. 
In order for the process to occur, masses of the matters which contribute to the anomaly 
diagrams must be smaller than the inflaton mass. Otherwise the matters decouple from 
the anomalies. For instance, when we consider the anomaly-induced decay into the SUSY 
breaking sector, since masses in the hidden quarks are expected to be of O (A), the decay 
takes place only for > A. 

So far, we have considered the anomalies of the SW-Kahler symmetry and sigma-model 
isometry. Since the process is an one-loop effect, there may be possible contributions from 
the counter term, depending on the underlying physics. Although we have assumed the 
conformal frame without the counter term at the cutoff scale in the above analyses, it 
can affect the decay rate, which is analogous to the anomaly-mediated SUSY breaking 



sccnano 



25|. 



Finally, let us comment on the inflaton decay into the SUSY breaking sector which 
involves the conformal dynamics. If the inflaton mass is above the scale of the violation 
of the conformal dynamics l , its decay into the SUSY breaking sector is expected to be 
suppressed. Actually, since the beta function vanishes above the scale, the decay induced 
by the SW anomaly is forbidden. At the same time, the contributions from the Kahler 
and sigma-model anomalies are implied to be suppressed m , once we notice that the SUSY 



breaking sector is sequestered from the other sectors by the conformal dynamics 28| . Then 
the inflaton field may not decay into the conformal SUSY breaking sector, and so, the models 
will be free from the gravitino production. 

A. Minimal Kahler Potential 

Let us explicitly show several examples of the spontaneous and anomaly-induced decays. 
The former decay depends on the Kahler potential of the inflaton and visible/hidden sectors. 
Let us first discuss the case of the minimal Kahler potential. We take the Kahler potential 



27] for a conformal theory of the SUSY breaking. 



See 

M. E. thanks K. -I. Izawa for discussions. 
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and the superpotential as 

K — (f>(f)* + QQ\ (29) 
W = W(<j>) + \MQQ + ^Y ljk Q l QiQ k , (30) 

where Q denotes the visible/hidden matters n . Then the total decay rate is the sum of the 
spontaneous and anomaly-induced decays, T = r( trce ) + r( anomal y)_ The former is given by 

r <-> „ t^W* U _ * + J^LM!, y . .W (311 

_ 167T M| Af£ * \ ml) 2567T 3 Mp ' u ' M p ( ' 

for fixed z, j and with i ^ j k. Here A^ 2 ) and A^ 3 ) denote a number of the final states. 
On the other hand, the anomaly-induced decay depends on the gauge structure. The decay 
rate is 

Wanomaly) _ NgO 2 2 (0) 2 m\ 

1 " 256tt 3 ^ G 1r) Ml Ml' {6I) 

In the above results, we have assumed that the inflaton mass is dominated by the SUSY- 
invariant mass term in the superpotential, and we have neglected the masses of the final 
states for the 3-body decay and the anomaly-induced decay. 



B. Sequestered Kahler Potential 



The next example is the Kahler potential with a sequestering form; 



K = -31oe 



1-\{H* + QQ* 



(33) 



with the superpotential (|30|) . Noting (T^) = (((f)) /3)5%, the rates of the spontaneous and 
anomaly-induced decays are 

144vr M p Mp \ ml J 

r (ano m al y) „ jV^O (ff m /o^N 

" 2304tt 3 Mp M P ' 1 ; 

where b is the beta function of the gauge symmetry, b = 3T G — T R . The tree-level decay 
arises via the mass term of Q, and the anomaly-induced decay is due to the SW anomaly, 



11 In addition, there may be a soft scalar mass in the Kahler potential. However it is irrelevant for the 
spontaneous and anomaly-induced decay processes. 
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while the spontaneous decay via the Yukawa coupling vanishes. We also find that the 
radiative effects associated to the Kahler and sigma-model anomalies cancel with each other, 
which is analogous to the cancellation of the AMSB effects to the gaugino mass 2^1 . 

The cancellation of the 3-body decay can be understood in the conformal frame. As 
was explained above, a part of the spontaneous decay is mediated by b m and the others are 
represented in terms of the chiral compensator field, $. The auxiliary field b m behaves as 
the gauge field of U(1)r of the super conformal symmetry. At the tree-level, this symmetry 
is preserved by assigning an U(1)r charge 2/3 for the chiral compensator. Then after a field 
redefinition, $Q — > Q, U(1)r charge vanishes for the operators which represent the 3-body 
final state of the decay, that is, 4 >% X''X k an d W^cj) 1 ^ (j)* k . Consequently, the decay mediated 
by b m cannot proceed via the Yukawa interaction. The other tree-level processes induced 
by the gravitational effects are also suppressed for the 3-body decay. They are obtained 
by multiplying $. Remembering that the anomaly- mediated SUSY breaking contributions 
to the soft trilinear couplings are absent at the tree level, $ does not contribute to the 
Yukawa interactions physically. Actually, $ becomes absent in the Yukawa interaction by 
the redefinition of the matte field, $Q — > Q- 

In addition to the gravitational effects discussed above, the inflaton decay may be induced 
by higher dimensional operators in the Kahler function, f2, which is defined as fl = —3e~ K / 3 . 
In fact, in contrast to the sequestered Kahler potential, there exist the higher dimensional 
terms in fl for the minimal Kahler potential. Thus the 3-body decays are allowed for the 
minimal case (see (13"T]) ). while they are absent in the sequestered one (see ( 1341) ) °. 



V. COSMOLOGY 



We now consider cosmological implications of the inflaton decay processes discussed in 
the preceding sections. One immediate consequence is that the reheating temperature Tr is 
bounded below; Tr cannot be arbitrarily low, since the inflaton decays into the visible sector 
through the top Yukawa coupling (See Eq. (fT8l) ). The other is the gravitino production 
from inflaton decay, which can occur through three different processes: (i) gravitino pair 

° In the Einstein frame, the cancellation can be seen explicitly by the field redefinition, e K / 6 Q — > Q, with 
K = K — {K) . This rescaling substantially corresponds to the transformation from the Einstein frame to 
the conformal one with respect to the interaction terms of the matters. 
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production; (ii) spontaneous decay at the tree level; (iii) anomaly-induced decay at the 
one-loop level. We will show how severely the gravitino production constrains the inflation 
models and SUSY breaking scenarios. 



A. Lower bound on the reheating temperature 

Let us begin with a relatively simple exercise. The supersymmetric SM sector contains 
the top Yukawa coupling in the superpotential as 

W = Y t TQH U} (35) 

where Y t is the top Yukawa coupling, and T, Q, and H u are the chiral supermultiplets of the 
right-handed top quark and left-handed quark doublet of the third generation, and up-type 
Higgs, respectively. In this section, we assume that the inflaton has the minimal Kahler 
potential for simplicity. The partial decay rate of the inflaton through the top Yukawa 
coupling is then 

Tt * T2^ |rt|2 ^ )2 ^' (36) 

where (<ft) and are VEV and the mass of the inflaton, respectively. The partial decay 
rate ( 1361) is derived from Eq. ( Tl8|) by noting ~ Y t 2 \ (0) | 2 and additional numerical factor 
6 coming from SU(3) x SU(2). The presence of the decay through the top Yukawa coupling 
sets a lower bound on the reheating temperature, T r. We define the reheating temperature 

as 

T « - (if) 4 yfc> ( 37 ) 

where g* counts the relativistic degrees of freedom, and denotes the total decay rate of 
the inflaton. Using Eqs. ( |36l) and ( 1371) . we obtain the lower bound on Tr, 



Tr > 1.9x10» GeV iKlf^V* (J$L 



'V (38) 



V 200 ) v 1015 GeV /v 1012 GeV 

Similarly the inflaton decays into the gluons and gluinos via the anomalies of SUGRA. 
One can estimate the decay rate from Eq. (128~1) as 

r S u(3) ^ (<P) 2 m% (39) 

where a s = gH^K denotes the SU(3)c gauge coupling constant. Substituting a s ~ 0.05, we 
can see r SU (3) is one order of magnitude smaller than T T . Therefore the spontaneous decay 
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into the visible sector is dominated by that through the top Yukawa coupling, unless the 
Kahler potential takes a specific form, i. e. the sequestered type (see Sec. UVBj) . 

We show the contours of the lower limit on Tr given by Eq. (I3"8j) in Fig. [21 together with 
typical values of {6) and for the single-field new 29|, multi-field new 30J, hybrid [3l| and 
smooth hybrid 32|, and chaotic 33j inflation models. We will discuss each inflation model 
later in this section. If the inflaton mass and the VEV (0) are too large, the reheating 
temperature may exceed the upper bound from the gravitinos produced by particle scattering 
in the thermal plasma. The cosmological constraints on the gravitino are summarized in 



Sec. IV Cl For more details, the reader should refer to Refs. 
gravitino, and Refs. 



34 



35 



37, 



361 ] for the unstable 



38l . 133 . |40| . l4ll . |42| for the stable one. For instance, the reheating 



temperature is necessarily higher than 10 6 GeV for the smooth hybrid inflation model, which 
is difficult to be reconciled with the gravitino of m 3 / 2 = 0(0.1 — 1) TeV [34J] and 10 eV ^ 



m 



3/2 



< 10 MeV U 



It is remarkable that the inflaton decays into the visible sector once it acquires a finite 
VEV; we do not need to introduce any interactions between the inflaton and the SM sector 
by hand in the Einstein frame p in order to induce the reheating. On the other hand, it 
may pose a cosmological problem at the same time. If the hidden sector also has a Yukawa 
coupling or includes the SW-Kahler/sigma-model anomalies, unwanted relics such as the 
gravitino may be directly produced by the inflaton. We will focus on the issue in the rest of 
this section. 



B. Gravitino Production 

We consider the gravitino production from the inflaton decay. To make our analysis 
simple and conservative, we assume that the inflaton has the minimal Kahler potential and 
does not have any direct couplings with the SUSY breaking sector in the superpotential. 
If we introduce possible couplings between the inflaton and the SUSY breaking field, the 
gravitino overproduction problem generically becomes severer. We also assume the DSB 
scenario with the dynamical scale A. Then the SUSY breaking field z usually has a scalar 
mass m z that is much larger than the gravitino mass. Although the precise value of m z is 

p Note that the interpretation of higher dimensional operators depends on a choice of the frame of SUGRA. 
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[GeV] 



+ : new(single);lTeV 
x : new(single);100TeV 
: new(multi) 

— : hybrid 

: smooth hyb. 

— : chaotic (w/o Z2) 



FIG. 2: Contours of the lower bound on Tr in units of GeV. We set g* = 228.75 and Yj = 0.6. 
For details of the models, see Sec. IV Ei 



model-dependent, it is expected to be of the order of A. Hereafter we simply assume q 

m z ~ A ~ y/m^. (40) 
We discuss the cases of > A and < A separately. 



1. The case oj ' m<f, < A 

As we have seen in Sec. Uli the inflaton decays into a pair of the gravitinos. The gravitino 
pair production is effective especially for a low-scale inflation model with < m z . The 
gravitino production rate is given by 



r 



(pair) 



(0>V 



(41) 



3/2 32tt 

for < m z ~ A. Here we have assumed the minimal Kahler potential with a soft scalar 
mass of z and (z) -C 1. The gravitino abundance is then 



3/2 



p(p air ) QT 

/3/2 O J- R 



q The scalar mass m z can be smaller than A. If this is the case, the pair-gravitino production will be 
affected. 
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~ 7 x io- 11 r^r 1 ( Tr y 1 ( w V r m * v m 

v2007 vl0 6 GeW ^KFGeV/ Uo 12 GeW ' 1 } 

It should be noted that the gravitino abundance is inversely proportional to Tr. This feature 
is to be contrasted to the thermally produced gravitinos, whose abundance is proportional 
to T R . 

2. The case of > A 

When the inflaton mass is larger than A, the gravitational effects discussed in Sec. IIHI 
and IIVI are important. If the SUSY breaking sector has Yukawa interactions, the inflaton 
decays into the sector via the operators. Besides, the anomalies of SUGRA induce the 
inflaton decay into the gauge boson and gauginos of the hidden gauge symmetries. Thus the 
hidden quarks and gauge bosons/gauginos are generally produced at the decay for > A. 

The hidden particles are energetic at the moment when they are produced. Since the 
reheating temperature Tr is bounded as Tr < A for almost entire region of the gravitino 
mass due to the thermal-gravitino production, the produced hidden particles do not reach 
thermal equilibrium. They instead form jets and hadronize by the strong gauge interactions, 
followed by cascade decays of the heavy hidden hadrons into lighter ones. The number of the 
hidden hadrons produced from each jet, which we call here as the multiplicity Nh, depends 
on the detailed structure of the hidden sector such as the gauge groups, the number of the 
matter multiplets, and a mass spectrum of the hidden hadrons. We expect N H to be in the 
range of 0(1 — 10 2 ). 

The hidden hadrons should eventually decay and release their energy into the visible 
sector, since otherwise they will easily overclose the universe. The gravitinos are likely to be 
produced in the decays of the hidden hadrons as well as in the cascade decay processes in 
jets. This happens, e. g. through the kinetic mbdngs of the hidden matters, and especially 
if z is a bound state of the hidden (s)quarks. Note that the goldstino is massless in the 
global SUSY limit and it is in the hidden sector with renormalizable couplings to other 
hidden (s)quark/gauge fields (and therefore hadrons). Thus, the goldstinos are expected to 
be produced by the hidden hadrons, though the precise production rate depends on details 
of the hidden sector. We denote the averaged number of the gravitinos produced per each 
jet as N 3 /2- Here we assume each hidden hadron produces one gravitino in the end, and use 
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the relation iV 3 /2 ~ r . 

The partial decay rates of the inflaton into the SUSY breaking sector are given by 



Eqs. 
tions 



43 



and 0321) . Although the DSB models do not always possess Yukawa interac- 
44 . 1451 ] . all the DSB scenarios necessarily involve the gauge interactions. From 
Eq. (I32p . the partial rate of the inflaton decay into the SUSY breaking sector is: 

Tdsb = -^£t(T^-T^r{<f>fm% (43) 

where the gauge coupling, the Dynkin indices, and the number of the generators are those 
of the hidden gauge symmetries. Multiplying the number of jets and A^/2, the gravitino 
abundance becomes 

Ys/2 = 2iV3/2 T74^' 

~ 9 x 10-^ (-^-Y 1 ( W V ( ^ V (44) 

" H 200/ vl0 6 GeW ^10 15 GeVj Uo 12 GeV; ' [ ) 

where we have defined £ = N 3 / 2 N^ hS) al{T G — T R ) 2 , which is roughly expected to be in the 
range of 0(l(r 2 ) to O(10). 

In the following numerical analysis, we take the anomaly-induced decay as a source of 
the gravitino production channel for > A. As one can see from Eqs. (I3"T|) and (I3"2"l) . the 
decay rate is roughly comparable to that of the spontaneous decay via the Yukawa coupling. 
Thus if one includes the tree-level decay into the analysis, the constraints become severer 
slightly, and the results in the followings do not change essentially. 



C. Cosmological Constraints on Gravitinos 



Before going further, here we briefly summarize the cosmological constraints on the grav- 
itinos, which will be used to put constraints on the inflation models later. 

There are tight constraints on the gravitino abundance from BBN if the gravitino is 



unstable 
38 



39 



34. 



35. 



36(| s , and from the dark matter (DM) abundance for the stable gravitino 37, 
The abundance of the gravitinos produced by thermal scatterings is 



In particular, if z is an elementary field and has a Yukawa coupling, the inflaton necessarily produces at 



least one goldstino by the decajLlhm 



For early works, see Refs. [46|, |47|, |48|, |49|, 15CJ, [51 



lgK, the, cca 



52, 
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related to Tr as [34j, [39] 



Y, 



(th) 



3/2 



1.9 x 1(T 12 
1 + 0.045 ln 




T 



10 10 GeV 
1 - 0.028 ln 



T 



R 



10 10 GeV 



(45) 



,10 10 GeV, 

where we have taken N = 3 for QCD and m g:j is the gluino running mass evaluated at 
T = Tr. Since the gravitino abundance Y^ is roughly proportional to Tr, Tr is bounded 
f r om above. 

Here we simply quote the bounds on Y 3 / 2 and Tr summarized in Ref. [7|. If the gravitino 
is light, it is likely the lightest SUSY particle (LSP) and therefore stable with the R-parity 
conservation. The bounds on Y 3 / 2 (and Tr) then come from the requirement that the 
gravitino abundance should not exceed the present DM abundance *: 



m 3/2 Y 3/2 < ^n DM < 4.4 x 10- 10 GeV, 



(46) 



where p c is the critical density, and we used Odm^ 2 ~ 0.12 at 95% CL. in the second 
inequality. The upper bound on T R can be obtained by substituting Eq. (I4"6"j) into Eq. (jl5|) 

as 



T R < 



O(100) GeV 
8 x 10 7 GeV 



for m 



3/2 



10' 



10 2 keV 



/ m 



m 



3/2 



for m 3/2 ~ 10- 4 - 10 2 GeV 



(47) 



V 300 GeV 7 VlGeV 

Note that we have conservatively neglected the contribution from the decay of the next- 
to-lightest SUSY particle. In the following analysis, we assume that the gravitino with 
a mass lighter than 10 2 GeV is the LSP and stable. When the gravitino is as light as 
m 3/2 ~ O(10) eV [55) ]. there are no constraints on T R} since the energy density of the 
gravitino would be too small even if the gravitino is thermalized. 



35|: 



On the other hand, if the gravitino is unstable, BBN puts severe constraints on Y 3 / 2 [34, 



Y 



3/2 



< 



1 x 10" 


-16 


- 6 x 10" 


-16 


for 


«13/2 


~ 0.1 


- 0.2 TeV 


4 x 10~ 


-17 


- 6 x 10" 


-16 


for 


"^3/2 


~ 0.2 


- 2 TeV 


7 x 10" 


-17 


- 2 x 10" 


-11 


for 


"^3/2 


~ 2 - 


10 TeV 


6 x 10- 


-13 


- 2 x 10" 


-12 


for 


"^3/2 


~ 10 - 


- 30 TeV 



{B h 



(4* 



' The gravitinos non-thermally produced by the inflaton decay can be a dominant component of DM, for 
certain values of the inflaton parameters 54 1 . 
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Y, 



3/2 



1 x 1(T 16 - 5 x 1(T 14 for m 3/2 
' < 2 x 1(T 14 - 5 x 10- 14 for m 3/2 
3 x KT 14 - 2 x 1(T 13 for m 3/2 
The corresponding upper bounds on T R are 

(1 - 4) x 10 6 GeV for m 3/2 - 

3 x 10 5 - 4 x 10 6 GeV for m 3/2 s 

5 x 10 5 - 1 x 10 8 GeV for m 3/2 s 

(3 - 10) x 10 9 GeV for m 3/2 - 



T, 



r ~ 



-0.1-1 TeV 
-1-3 TeV 
-3-10 TeV 

0.1-0.2 TeV 
0.2 - 2 TeV 
2-10 TeV 
10 - 30 TeV 



(B h ~ 10- 3 ). (49) 



(B h ~ 1), 



(50) 



Tr < 



10" 



(51) 



1 x 10 6 - 3 x 10 8 GeV for m 3/2 ~ 0.1 — 1 TeV 
(1 - 3) x 10 8 GeV for m 3/2 ~ 1 - 3 TeV (B h 
^ 2 x 10 8 - 1 x 10 9 GeV for m 3/2 -3-10 TeV 

For the heavy gravitino of mass ^ 30(10) TeV, no stringent constraints are obtained from 
BBN. However, another constraint comes from the abundance of the LSP produced by the 
gravitino decay. Since the gravitino life time is rather long, the produced LSPs will not 
annihilate with each other. Thus the upper bounds on Y 3 / 2 and Tr read 

m LSP V 3/2 < 4.4 x 10- 10 GeV, (52) 



and 



,io ( "^lsp 
H ~ - ' JXi ° vI^GeV 



Tr < 2.5 x 10 1 



GeV, 



(53) 



where ?71lsp denotes the mass of the LSP. 

As is well-known, all the above constraints have been usually applied for the gravitinos 
from the thermal production. Since the gravitinos are also non-thermally produced in infla- 
ton decay, we obtain further constraints on T R , m 3 / 2 , (0) and by requiring the abundance 
of the non-thermally produced gravitinos (j42]) and (jHj) to satisfy (SED , P*)) . (H9|) . or (JS2J) - As 
we will see later, these new constraints drive (some part of) the high-scale inflation models 
and the gravity mediation into a corner. 



D. Constraints on Inflation Models and SUSY breaking 



Now we would like to derive constraints on the inflation and SUSY breaking models, using 
the non-thermal production of the gravitinos discussed above together with the thermal 

23 



process. 

In Fig. El we show the constraints on the inflaton mass and VEV for m 3 / 2 = 1 GeV, 1 TeV, 
and 100 TeV, together with typical values of the inflation models. We discuss each model in 
the next subsection. The region above each solid line is excluded. We find that in the case of 
m 3/2 = 1 TeV with = 1, all the inflation models shown in the figure are excluded. For the 
gravitino mass lighter or heavier than the weak scale, the constraints become relaxed. The 
inflaton mass and its VEV depend on the inflation models. Generically speaking, for larger 
and ((p), the constraints become severer, simply because more gravitinos are produced 
by the inflaton decay (see (j4*2l) and (PH]) ). On the other hand, if the inflaton is charged 
under some symmetries, its VEV becomes suppressed or even forbidden especially when 
the symmetry is exact at the vacuum. Then the bounds can be avoided for such inflation 
models. This is the case of the chaotic inflation model with a discrete symmetry (note that 
the chaotic inflation model shown in Fig. Elis that without such a symmetry). 

The solid lines which denote the constraint are jaggy at an intermediate value of m^. 
This is because the dominant production channels of the gravitinos changes. In the right 
side, the gravitinos are produced by the spontaneous and anomaly-induced decays, while 
the inflaton directly decays into a pair of the gravitinos in the left side. Note that we have 
assumed (flO]) and £ = 1 for simplicity. 

In Fig. El we have set Tr to be the highest value allowed by the constraints. As mentioned 
before, the abundance of the non-thermally produced gravitinos is inversely proportional to 
Tr, which is different from that of the thermally produced one (see f T4"5l) ). If Tr takes a 
smaller value, the constraints becomes severer. Thus, the bounds shown in Fig. El are the 
most conservative ones. Note that one may have to introduce couplings of the inflaton with 
the SM particles to realize the highest allowed reheating temperature. 

Instead, taking Tr as a free parameter, we show the constraints on the — Tr plane, 
for m 3/2 = 1 GeV, 1 TeV, and 100 TeV with a fixed (0) = 10 15 GeV in Fig. PJ The reheating 
temperature is bounded from above due to the thermal production of the gravitino. It is 
remarkable that we have lower bounds on Tr due to the non-thermal processes. In the 
figure, we incorporated the spontaneous decay via the top Yukawa interaction, which also 
provides a lower bound on Tr (see (138|) ). One can see that the lower bound on Tr becomes 
severer for larger m^. 

In Fig. El we show constraints on the plane for several values of (<p) = 10 12 , 10 15 , 
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+ : new(single);lTeV 
x : new(single);100TeV 
: new(multi) 

— : hybrid 

: smooth hyb. 

— : chaotic (w/o Z2) 



[GeV] 

FIG. 3: Constraints from the gravitino production by the inflaton decay, for m 3 / 2 = ITeV 
with Bh = 1 (caseA), m 3 / 2 = 1 TeV with B^ = 10~ 3 (caseB), 7713/2 = lOOTeV (caseC), and 
m 3/2 = 1 GeV (caseD). The region above the solid (gray) line is excluded for each case. For 
rricf, > A, we have used the anomaly-induced inflaton decay into the hidden gauge/gauginos to 
estimate the gravitino abundance, while the gravitino pair production has been used for < A. 
Since Tr is set to be the highest allowed value, the constraints shown in this figure are the most 
conservative ones. 



and 10 18 GeV. The dashed (pink) line represents = A. For the inflaton mass above 
the dashed (pink) line, the spontaneous and anomaly-induced decays of the inflaton produce 
the gravitinos, while the pair production is dominant below the dashed (pink) line. We have 
set Tr to be the highest value allowed by the constraints as we did in Fig. [3j We find that 
the inflaton mass cannot be too large, especially for m 3 / 2 around the weak scale. It is also 
noticed that the constraint becomes severer as (0) increases, since the upper bound on 
is proportional to for fixed m 3 / 2 . 

Finaliy let us illustrate how much the problem becomes severer in the case of the SUSY 
breaking models with an elementary singlet z. Such a field is needed to give sizable masses 
to the gauginos in the simple version of the gravity-mediation Q, Q (see also footnote [5J) . 
For an inflaton with < A, the pair gravitino production occurs as described before. In 
particular, since z is singlet at the cutoff scale, there is a priori no reason to forbid such 
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FIG. 4: Constraints from the gravitino production by the inflaton decay, for 1713/2 = 1 GeV (left- 
upper), m 3 / 2 = ITeV with Bh = 1 and 1CT 3 (right-upper) , m 3/ / 2 = 100 TeV (bottom). We have 
set {(j)) = 10 15 GeV. The region surrounded by the solid line is allowed for each case. 



an interaction as 5 K ~ \(p\ 2 (z + z*). Then there generically exists a large kinetic mixing 
with the inflaton, and so, the gravitino production rate becomes too large, which is given 
by Eq. with c ~ 1. With such a large gravitino production rate, most of the inflation 
models with < A are excluded, e. g. unless the inflaton VEV vanishes due to some 
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FIG. 5: Constraints from the gravitino production by the inflaton decay, for (</>) = 10 12 , 10 15 and 
10 18 GeV. The region above the thick solid line is excluded. We also show the constraint for the 
unstable gravitino with Bh = as the thin (blue) line. For the region above the dashed (pink) 
line, we adopt (|44j) . while (I42j) is used for the region below the dashed line. Since Tr is set to be 
the highest allowed value, the constraints shown in this figure are the most conservative ones. 



symmetries. Even for > A, the pair gravitino production occurs effectively. In fact, one 
expects that 5K ~ \(p\ 2 zz/2 + h. c. generally exists. The gravitino pair production is then 
given by (jSJ) with c ~ 1. Thus the gravitino abundance increases by O(10 2 ) compared to 
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that from the spontaneous and anomaly-induced decays with £ = 1 for > A (see (1421) and 
(j44j) ) u . Besides, such a z field may be displaced away from its potential minimum during 
the inflation, forcing the cosmological scenario to be more problematic (see the footnote [cj) . 
Thus the SUSY breaking models with the elementary singlet z are strongly disfavored from 
the cosmological points of view. 



E. Inflation Models 



In this subsection we give a brief review on the representative inflation models plotted in 
Figs. [2] and [3j For details on the models, the readers should refer to the original literatures. 



1. Single-field inflation model 



As a concrete example, here we study the new inflation model [29l . l56l . 1571 ] . In the new 
inflation model, the Kahler potential and superpotential of the inflaton sector are written 
as v 

Kfatf) = |0| 2 + ^|0| 4 , 

W W = v ^--J- l ^\ (54) 

where the observed density ffuctuations are explained for v — 4 x 10~ 7 (0.1 / g) 1 ^ 2 and k < 
0.03 in the case of n = 4 After inflation, the inflaton takes the expectation value 
(0) ~ (v 2 / g) l l n . In this model the inflaton mass is given by ~ nv 2 j (0), and the gravitino 
mass is related to v as m 3 / 2 — nv 2 (0) /(n+ 1), since the inflaton induces the spontaneous 
breaking of the i?-symmetry. 

In the case of n = 4, the inflaton parameters are ~ 4 x 10 9 GeV and (0) ~3x 10 15 
GeV for m 3/2 = 1 TeV, while m^~2x 10 10 GeV and (0) ~ 1 x 10 16 GeV for m 3/2 = 100 TeV. 
Note that m 3 /2 <C 1 TeV cannot be realized unless g ^> 1. From Fig. [31 we can see that the 
new inflation model is excluded for m 3 / 2 = 1 TeV with B^ — l, while it is below the bound 
for m 3/2 = 100 TeV. 



u The spontaneous decay at the tree level and the anomaly-induced one into the SUSY breaking sector are 

not much affected by the presence of such an elementary singlet z. 
v The gravitino abundance in the text remains virtually unehanged in the presence of the quartic coupling 

in the Kahler potential. 
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2. Multiple-field inflation model 



Next we consider an inflation model with multiple fields. Among many multiple-field 
inflation models proposed so far, there is an important class of models described by the 
following superpotential: 

W(M) = <f>m, (55) 

where f(if>) is a function of ip. The potential minimum in the global SUSY limit is located 
at 

(0) = 0, 

m = ^ (56) 

where tpo satisfies f(ipo) = 0. Note that the true minimum is slightly displaced from (1561) . 
once the SUSY breaking field is taken into account H]- 

For instance, the above class of the models includes a new inflation model [3^] and a 
hybrid inflation model [311, described by 



W(M) = 0(^-^1, (57) 



where /i determines the inflation energy scale and M is an effective cut-off scale. In the new 
inflation model ip plays a role of the inflaton, while is the inflaton in the hybrid inflation 
model. 

The inflaton fields <fi and ip have almost the same masses, 



m ~ 



e G / 2 V^ , (58) 



which are assumed to be much larger than the gravitino mass. It should be noted that 
(p and ^ (and/or %p) almost maximally mix with each other to form the mass eigenstates 
due to the almost degenerate masses [7]. To see this one should note that the difference 
between the diagonal components of the mass matrbc is small: \M^ — M^\ = 0(m|/ 2 ), while 
the off-diagonal component is relatively large: ML, = 0(m3/ 2 m^), resulting in the almost 
maximal mixing between and if}'. Similar mixing may occur between and ip. This 
mixing is effective at the inflaton decay, since the Hubble parameter at the decay should be 
smaller than 0(m 3 / 2 ) to satisfy the bounds from the thermally produced gravitinos. The 
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mass eigenstates are obtained after taking account of the (almost) maximal mixing between 
4> and 

<P± ^ ^^p- (59) 
The mass-eigenstates have the mass given by (j58p and the effective VEV (f±) given by 
ipo/\/2 unless there is cancellation. 



2- A. New inflation model 



The new inflation discussed above is also realized for 30j 



k = i0r+i^r+^i0i 4 +A; 2 |0i 2 i^r+|i^i 4 , 

W = (f)(v 2 -g^), (60) 

in which the inflaton is ip, while <fi stays at the origin during and after inflation. If one defines 
k = k,2 — 1, the scalar potential for the inflaton i/j becomes the same as the single-field new 
inflation model, although the gravitino mass is not related to the inflaton parameters. After 
the inflation ends, the energy of the universe is dominated by the oscillation energy of i/j. 
Although 0o is suppressed compared to ipo, the effective VEV is given by tpo/\/2, since 
and ip^ almost maximally mbces with each other in the vacuum. Thus the constraint on this 
model is comparable to that on the single-field new inflation. We plot the values of and 
(<p±) for g = 10~ 4 — 1 and k = 10~ 4 — 10 -1 ' 5 with the e-folding number = 50 in Figs. [2] 
and [3l The (multi-field) new inflation model is excluded for m 3 / 2 = 1 TeV with — 1, 
while it is allowed for m 3 / 2 = 1 GeV and 100 TeV. 

2-B. Hybrid and Smooth hybrid inflation models 
The hybrid inflation model contains two kinds of superfields: one is which plays a 



role of inflaton and the others are waterfall fields ip and ip 31]. After inflation ends, <fi as 
well as oscillates around the potential minimum and dominates the universe until the 
reheating. 

The superpotential W(<f>,if},i^) for the inflaton sector is 

iy(0,^,^) = 0(/i 2 -A^), (6i) 
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where ip and ip are assumed to be charged under U(l) gauge symmetry. Here A is a coupling 
constant and /i is the inflation energy scale. The potential minimum is located at (0) = 
and (?/>) = (?/>) = /i/ \f\ in the SUSY limit. For a successful inflation, /i and A are related as 
/i~2x lO^A 1 / 2 for A > 1(T 3 , and /i ~ 2 x 10~ 2 A 5 / 6 for A < 10~ 3 . 

Due to the D-term potential one linear combination of ip and ip, given by ip^ = — 
tp)/y2, has a large mass of ~ g (ip) (g denotes the gauge coupling), while the other, = 
(i/i + $)/V2 has a mass equal to that of 0: m^+) = = \/2X(ip). It is the latter that 
(almost) maximally mixes with <\> to form mass eigenstates. Note that VEV of is equal 
to y/2(tf>). 

For A ~ 10" 1 - 10^ 5 

and (<p±) = fi/y/X ~ O(10 15 )GeV . From Fig. El one can see the hybrid inflation model is 
excluded by the gravitino overpro duet ion for m 3 / 2 = 1 TeV with B 3 / 2 — 1- For m 3 / 2 = 1 GeV 
and 100 TeV, the constraints become slightly mild, but a certain fraetion of the parameter 
space is still excluded. The allowed parameter space corresponds to A ^ 10~ 2 . Note that 



5^ we obtain /i~8x 10~ 4 - 1 x 10" 6 , ~ 10 15 - 10 10 GeV, 



eads to almost scale-invariant 
2|. It is possible to make the 



the parameter space allowed by the gravitino produetion 
power spectrum, which is disfavored by the WMAP data 
scalar spectral index ru smaller than 1 by introducing non-renormalizable interaetions in 
the Kahler potential 59l. l60|. 

Here we comment on interesting observation concerning the spectral index and the cosmic 
string. In this type of hybrid inflation, cosmic strings are formed after inflation because ip 
and ip have U(l) gauge charges. As is well known, the cosmic strings contribute to the 
density fluetuations. Including the effects of the cosmic st ring makes the spectral index 
n s between 0.98 and 1 compatible with the WMAP data 6jJ, 162] , if the tension of the 
cosmic string is Gfi = O(10 -7 ). According to Ref. 6jJ, this corresponds to the region with 
A ~ 0(1O~ 3 - 10~ 2 ). Interestingly enough, the region is just below the constraints from 
the gravitino produetion in the case of = 1 GeV and 100 TeV w . This means that, 
for that region, the gravitino non-thermally produced by the inflaton decay may account 
for the present DM abundance 54j, if the gravitino is stable. For the unstable gravitino 



of a mass m 3 / 2 ^ O(10) TeV, Wino-like LSP produced by the gravitino decay may be the 



Including the soft terms, the inflaton dynamics is somewhat modified, and correspondingly the inflaton 
parameters are slightly changed, especially if the gravitino mass is heavy 63 1. 
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dominant component of DM. Moreover, since the required tension of the cosmic string is 
relatively large and is close to the present observational upper bound, one may be able to 
discover the cosmic string in the future observations. Since the inflaton mass and the VEV 
are small, it is difficult to realize the non-thermal leptogenesis via the spontaneous decay 
(see Eq. (1151) ). However, one can naturally incorporate the non-thermal leptogenesis into 



the hybrid inflation model by identifying the U(l) symmetry with a U(1)b-l symmetry [64]. 

Next let us consider a smooth hybrid inflation model j^j], which predicts the scalar 
spectral index as n s ~ 0.97, which is slightly smaller than the simple hybrid inflation model. 
The superpotential of the inflaton sector is 

nM,t)=4>U-^)- (62) 



The VEVs of ip and ip are given by (ip) = (ip) = (/iM"" 1 ) 1 '' 1 , and we assume that ip = ip 
always holds due to the additional D-term potential. Then one of the combination, ip( + ' = 
(ip + ip)/ y/2, almost maximally mixes with <p> to form the mass eigenstate of a mass = 
y/2nn 2 / {ip). For n = 2 we obtain /i ~ 4x i(T 4 -9 x 1(T 5 , and m^~lx 10 14 -6 x 10 14 GeV. 
From Fig. [31 one can see that the smooth hybrid inflation model is excluded for a broad 
range of m 3 / 2 . 



2- C. Chaotic inflation model 

A chaotic inflation is realized in SUGRA, based on a Nambu-Goldstone-like 

shift symmetry of the inflaton chiral multiplet <p 33j. Namely, we assume that the Kahler 
potential K(cp, (p') is invariant under the shift of (p, 

<p^<p + iA, (63) 

where A is a dimensionless real parameter. Tirus, the Kahler potential is a function of (p + (p^; 
K (<p, (p') = K (<p + (p') = c (0 + (p') + \((p + (p^) 2 + ■ • •, where c is a real constant and must 
be smaller than 0(1) for a successful inflation. As opposed to the other inflation models, 
this model allows a linear term in the Kahler potential. The coeflicient c corresponds to the 
inflaton VEV in the other models. If there is no other symmetry such as a Z 2 symmetry, 
there is no reason to expect that c is much smaller than unity. 
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We identify the imaginary part of (p with the inflaton field ip = y/2 Im [<p]. Moreover, we 
introduce a small breaking term of the shift symmetry in the superpotential in order for the 
inflaton (f to have a potential: 

W((f>,ip) =m<j)ip } (64) 

where we introduced a new chiral multiplet ip, and m ~ 2 x 10 13 GeV determines the inflaton 
mass. 

One might suspect that it is only the real component of <fi that can decay into the 
gravitinos, since the shift symmetry dictates that the only real component (<f) + <fi*) appears 
in the Kahler potential. However, it is not surprising that this is not the case, since the 
decay amplitude is proportional to powers of the large SUSY mass m that explicitly violates 
the shift symmetry. 

We plot the chaotic inflation model with c = 0.1 — 1 in Figs. [2] and [31 One can see that 
it is excluded for almost entire values of m^n (except for m 3 / 2 % O(10) eV). Note however 
that one can avoid the constraints by assuming an approximate Z2 symmetry to suppress c. 



F. Possible solutions to the gravitino problem 



Here let us briefly mention possible solutions to the gravitino overpro duet ion problem. As 
mentioned above, one solution is to postulate a symmetry of the inflaton. If the symmetry 
is unbroken at the vacuum (or if the breaking of this symmetry is small), the VEV of the 
inflaton, (0), is zero (or suppressed). As the gravitino produetion rate is proportional to 
(4>) 2 , one can avoid the gravitino overpro duet ion for such inflation models. Note however 
that, if the symmetry is exact, the visible matter fields as well must be charged under 
the same symmetry, since otherwise the inflaton cannot decay into the visible seetor. This 
solution can be achieved e. g. in the chaotic inflation model; one can assign Z2 symmetry on 
the inflaton 20]. Also, t here are inflation models in which the inflaton is identified with the 



MSSM fields 
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61 



691 ] or the right-handed sneutrino [70]. By similar reasoning, the 



gravitino overpro duet ion from the inflaton decay can be avoided in these models. 

So far, we have assumed that there is no late-time entropy produetion after the inflaton 



decay. If huge entropy produetion 







721 ] oecurs after the reheating of the inflaton, any pre- 



existing gravitinos are diluted. However, since it also dilutes the pre-existing baryon number, 



one may have to generate the baryon asymmetry after the entropy produetion 



73 



75 



76] 
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Note also that the decay processes discussed in this paper can be applied to any scalar fields, 
and so, the scalar field that induces the large entropy may produce the gravitinos again. One 
has to make sure that this does not happen. 

Another solution is to assume that the gravitino mass is either extremely heavy or ex- 
tremely light. If the gravitino mass exceeds O(10 6 ) GeV, the gravitino may decay before 
the decoupling of the LSP particle. To realize this, however, one has to contrive a set-up 
in which anomaly-mediation is suppressed. If we stick to the gravitino mass smaller than 
lOOTeV on the basis of naturalness, the gravitino problem sets severe bounds on the high- 
scale inflation models as seen above. On the other hand, if the gravitino mass is lighter 
than 0(10) eV, the gravitino is thermalized and its cosmological abundance is negligibly 



small 



55j. So, the gravitino problem is absent for such a very light gravitino. Note that, 



in order to ameliorate the gravitino problem instead of solving it completely, one does not 
have to go to such extremes. For the gravitino mass moderately lighter or heavier than the 
weak scale, the gravitino problem is relaxed especially for the low-scale inflation models (see 

Fig.E]). 

If the inflaton mass is quite light, one can evade the constraints as one can see from Fig. [31 
However, it should be noted that the reheating temperature is set to be the highest allowed 
value in Fig. [3j For a lower reheating temperature, the constraints become severer. This 
means that, one may have to introduce relatively strong couplings of the inflaton to the SM 
particles in order to realize the reheating temperature adopted in Fig. [3l For instance, if the 
inflaton has only the dimension-five couplings to the visible sector suppressed by the Planck 

3 /2 

scale, the reheating temperature Tr is proportional to ml , which makes the gravitino 
abundances (j4"2j) and (1441) rather insensitive to the inflaton mass. Then, one cannot evade 
the gravitino problem simply by changing the inflaton mass x . 

Lastly, let us comment on the SUSY breaking sector which involves the conformal dy- 
namics. As mentioned in Sec. IIV| the inflaton field with m$ ^> A may not decay into 
the conformal SUSY breaking sector. Then the gravitino production is suppressed. This 
solution is appealing because one does not have to impose non-trivial constraints on the 
inflation models or on the thermal history of the universe. The only requisite is the con- 
formal dynamics in the SUSY breaking sector, which has its own phenomenological virtues 



x This is one of the reasons why the moduli-induced gravitino problem (4] is quite difficult to be solved. 
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independently of the gravitino overproduction problem. Furthermore, it may naturally lead 
to the successful non-thermal leptogenesis scenario, which will be discussed elsewhere [77 ]. 



G. Comments on preheating 



Here we would like to mention the effects of the preheating (78j, i. e., the non-perturbative 
inflaton decay process, which we have not taken into account so far. The non-thermal 
gravitino production from the inflaton decay has somewhat checkered history; it was once 
claimed that the gravitinos were non-thermally produced during preheating [791 ]. but it was 
later concluded that the inflatino, instead of the gravitino in the low energy, was actually 
created [80 ] . 



Since the inflatino decays much earlier than the BBN epoch |8l|], the non- 
thermal 'gravitino' (actually, inflatino) production turned out to be harmless. However, 
as discussed so far, we have found that the gravitinos are generically produced by the 
perturbative decay processes. 

Another concern is whether our results are modified by including the effects of the pre- 



heating 82]. We believe that our arguments on the gravitino overproduction problem are 
robust and they are not essentially modified even if the preheating occurs. First of all, we 
would like to emphasize that the decay processes discussed so far are perturbative ones, 
and therefore they are always present. On the other hand, it crucially depends both on the 
global structure of the inflaton potential and on the couplings of the inflaton to matter fields 
whether the preheating occurs and how efficiently it proceeds. 

Let us assume that the preheating actually occurs and it proceeds quite efficiently with- 
out any back reaction, i. e., the inflaton transfers most of its energy into other particles soon 
after the inflaton starts oscillating. This corresponds to the instantaneous reheating, which 
generically leads to the overproduction of the gravitinos due to particle scatterings, instead 
of the non-thermal production. The latter is suppressed in this case since the reheating 
temperature will become high (see (|42p and fj44l) ). Therefore, for the most inflation mod- 
els, such an efficient preheating should not occur, since otherwise too many gravitinos are 
produced by the conventional thermal scatterings. 

On the other hand, if the preheating is not so efficient, then one has to take account 
of the back reaction processes, and the preheating typically ends at a certain point. The 
reheating of the universe is induced by the perturbative decay of the inflaton in the end, 
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and the gravitinos are generically produced by the decay. Thus it is unlikely that the 
preheating solves or ameliorates the gravitino problem; one has to contrive a model in 
which the preheating proceeds quite efhciently, but the gravitinos are not produced by the 
scatterings of the decay products. Specifically, the decay products should not reach thermal 
equilibrium. 

VI. CONCLUSION 

The present observational data on the CMB and the large-scale structure, together with 
the strong theoretical motivation to resolve severe problems in the Standard big bang cos- 
mology, have led us to believe that the universe underwent an inflationary epoch at an early 
stage. While there are many inflation models (called as "the inflationary zoo"), we still do 
not know which inflation model is realized in nature. The study on the density fluctuations 
such as isocurvature perturbations, non-gaussianity, tensor-mode, and their effects on the 
CMB power spectrum is quite useful, but is not enough at present to pin down the inflation 
model. This is partly because of our ignorance of thermal history of the universe beyond 
the Standard big bang theory, e. g., how the inflaton reheats the universe. 

In this paper, we have investigated the inflaton decay processes in the supergravity. In 
particular, we have shown that the gravitinos are generically produced in the inflaton decay. 
There are three different processes for the production. One is the direct production of a pair 
of the gravitinos. This is effective especially for low-scale inflation models. The other two 
are due to the inflaton decay into the SUSY breaking sector; the spontaneous decay at the 
tree level and the anomaly-induced one at the one-loop level. Those non-thermally produced 
gravitinos set tight constraints on the inflation models, together with the constraints from 
thermally produced gravitinos. Indeed, these two constraints are complementary in a sense 
that the dependence on the reheating temperature is different. For higher Tr, more grav- 
itinos are thermally produced, while the non-thermally production becomes important for 
lower T R . They also depend on the SUSY breaking scenarios. In fact, almost all parameter 
space for the inflaton is excluded for the gravity-mediated SUSY breaking scenarios. 

Apart from the gravitino productions, the inflaton naturally decays into the visible sector 
especially through the top Yukawa coupling and SU (3) c gauge interaction. Thus one does 
not need to introduce any direct ad hoc couplings by hand in order to induce the reheating. 
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The above studies may provide us with a breakthrough toward the full understanding of 
the inflationary universe. In addition to the Standard analysis on the density fluctuations, 
the inflation models in supergravity are subject to the constraints due to the (non)-thermally 
produced gravitinos. Whether a consistent thermal history after inflation is realized now 
becomes a new guideline to sort out the inflationary zoo, and hopefully it will pin down the 
true model, together with data in the future collider experiments such as LHC. 
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